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Abstract 
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1 Introduction 

For a,/3 > —1, the Hahn polynomials are explicitly given by 



Qnix;a,/3,N) := 3F2{-n, -x,n + a + l3 + 1; -N, a + 

k{-x)k{n + a + /3 
{-N)kia + l)kk\ 



E {-n)k{-x)k{n + g + /3 + l)fc _ 



k=0 

see [2, p. 174] and [6]. These polynomials are orthogonal on the discrete set {0, 1, • • • , A^} with respect to 
the weight function 

p[x;a,/3,N) = (^jsf _ ' x = 0, 1, 2, • • • , iV. (1.2) 

More precisely, we have 

N 

Qn{k; a, (3, N)Qm{k; a, (3, N)p{k; a, /?, N) = /iAr,„5„,m, n, m = 0, 1, • • • , A^, (1.3) 

k=0 

where 

, ^ (-l)"(n + a + /3 + l)jv+i(/3 + l)nn! 
"^'^ {2n + a + f3 + l)(a + l)„(-iV)„iV! ' ^ ' ' 

see [9, p. 462]. 
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Formula (1.3) tells us that these polynomials are orthogonal on an unbounded interval as n — )• oo. We 
now introduce a rescaling so that these polynomials become orthogonal on a bounded interval. Let X^v 
denote the set defined by 

k -\- 1/2 

Xn ■■= {xN,k]k=o^ where xj^^k ■= — — • (1-5) 
The X7v,fc's are called nodes and they all lie in the interval (0, 1). Also, let 

PnA^) ■= QniNz - 1/2; a,P,N- 1) (1.6) 

and 

w{z) := N-'^-^a\p\p{Nz - 1/2; a,P,N - 1). (1.7) 

It is readily seen that the polynomials PN,n{z) are orthogonal on the nodes x^^j. with respect to the weight 
WN,k ■= 'w{xN,k); that is, 

N-l 

PN,n{xN,k)PN,m{xN,k)wN,k = h*^,n^n,mi (1.8) 

fc=0 

where h*^^^ = A^-°-^a!/3!/i^_i,„. 

For further properties and applications of the Hahn polynomials, we refer to Beals and Wong [2] and 
Karlin and Mcgregor [6]. In 1989, Sharapudinov [13] studied the asymptotic behavior of Qnix; a, (3, N), 
when the degree n is substantially smaller than N. His main result is an asymptotic formula for 
Qnix;a,p,N) with n = 0{N^/'^) as n — )■ oo; which involves the Jacobi polynomial pji^'°'\t), where t 
is related to x via the formula x = {N — 1)(1 + 1)/2. Recently, Baik et al. [1] used the Riemann-Hilbert 
approach to investigate the asymptotics of discrete orthogonal polynomials with respect to a general 
weight function. Their results are very general, but it is difficult to use them to write out explicit for- 
mulas for specific polynomials. With regard to the Hahn polynomials, they only considered the case of 
varying parameters; that is, a = NA and /3 = NB, where A and B are fixed positive numbers. Moreover, 
their results are more local in nature; that is, one needs more asymptotic formulas to describe the behavior 
of these polynomials in the complex plane. 

The purpose of this paper is to study uniform asymptotic behavior of the Hahn polynomials as n — )• oo, 
when the parameters a and f3 are fixed and the ratio n/N is a constant c € (0, 1). Our approach is based 
on a modified version of the Riemann-Hilbert method introduced by Deift and Zhou [4]. This version has 
been used successfully to obtain globally uniform asymptotic expansions of several discrete orthogonal 
polynomials; see [3], [10] and [15]. The presentation of this paper is arranged as follows: In Section 2, 
we review some preliminaries, including weak asymptotics of the zero distribution and the formulation 
of the basic interpolation problem Y(z). In Section 3, we transform this interpolation problem into an 
equivalent Riemann-Hilbert problem (RHP). In Section 4, we give the matrix transformation T{z) that 
normalizes the RHP for R{z) presented in Section 3 by using a function g{z), which is related to the 
logarithmic potential (g- function) of the equilibrium measure. The auxiliary functions (j){z) and (j)*{z) are 
also studied in Section 4. In Section 5, we factorize the jump matrix in the RHP for T{z). The final 
results are presented in Sections 6 and 7, where we divide the whole complex plane into three regions 
and in each region we derive a uniformly asymptotic formula. As a special case of our result, we give in 
Section 8 the asymptotic formula of Qnix; a, /3, N) for fixed values of x. 
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2 Preliminaries 



2.1 The equilibrium measure 

It is known that the zero distribution of discrete orthogonal polynomials is related to a constrained 
equilibrium problem for logarithmic potential with an external field '■p{x)^ which is defined by the formula 

ip{x) := V {x) + I \n\x -y\p^{y)dy, a; €(01,02), (2.1) 

J ai 

where V{x) and p^{x) are real-analytic functions in a neighborhood of the closed interval [oi, 02], and p^{x) 
is also a density function of the nodes; see [3] and [1, p. 26]. In our case, V{x) = 1— xlogx— (1— x) log(l— x), 
p^{x) = 1 and [ai, 02] = [0, 1]. To see how the function V{x) is derived, we recall the nodes XN,n and the 
weights wiq^n = w{xN^n) given in (1.5) and (1.8), respectively, and write 



7V-1 

WN,n = 



-Q (2.2) 



m=0 



Prom (2.2), an explicit formula can be given for VN{x]^^n)- Indeed, since x^^n = {n + \)/N and 



N-l ^r„i 



JJ \xN,n — XN,m\ 



n\iN -n-l)V 

m=0 ^ ' 

we have 

1 1 N^- 

yN{xN,n) = -TT \.Og w{xN,n) + — log ■ 



N ^ ' ' ' N ^ n\{N -n - l)V 

Since n = Nxj^^n ~ ^1 by replacing XK[,n and n by x and Nx — ^, respectively, in the above equation, we 
obtain 

1 1 N^^^ 

Vn{x) = ——logwix) + —\og T- T— (2.3) 

^' N ^ N {Nx - ^)\{N - Nx - ^)\ ^ ^ 

for X G (0, 1). Using Stirling's formula, it can be shown that 

1 iv^-i / 1 

TT log -, TTT, TTT = 1 - X log X — (1 — x) log(l — x) + O [ — 

N {Nx -ly.{N - Nx B \ J &\ ) 

as — 7- oo. From (1.7), we also have 

,T{Nx + a + \) T{N{l-x) + 13+\) 



w{x) = iV^°^^- 



T{Nx + \) T{N{l-x) + \) 
~x°(l-x)'^ asiV^cx). 

Thus, it follows from (2.3) that 



Vn{x) = 1 -xlogx- (1 - x)log(l -x) + 0[^] (2.4) 
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as — )• oo. This approximation formula suggests the definition of the function 



V{x) = 1 — X log X — {I — x) log(l — x), (2.5) 

and we have 

VNix) = V{x) + o(^-^y (2.6) 

Observe that in our case, we also have <f{x) = 0. 

The measure fj,{x)dx that minimizes the following energy functional 

Ec[l-i] = —c / / In \x — y\fi{x)dxfi{y)dy + / ip{x)fi{x)dx, (2-7) 
Jo Jo Jo 

where fi{x) satisfies the upper and lower constraints 

< fi{x) < 1/c (2.8) 

and the normalization condition 

1 

H{x)dx = 1, (2.9) 



is called the equilibrium measure. It follows from (2.8) that there is an upper bound for ^(x), which 
does not appear in the continuous case. This fact is a crucial difference between discrete orthogonal 
polynomials and continuous orthogonal polynomials. The equilibrium measure fi{x)dx divides the interval 
of orthogonality into subintervals of three possible types: (1) achieving the lower constraint; (2) attaining 
the upper constraint; (3) not reaching the upper or the lower constraints. We call the open intervals of 
type (l)-(3), Voids, Saturated Regions and Bands, respectively. For reference, we refer to [1, p. 19]. 

In the existing literature, the equilibrium measure is usually obtained by solving a minimization 
problem of a certain quadratic functional (cf. [1], [3]). Here, we prefer to use the method introduced by 
Kuijlaars and Van Assche [7]. 

Following [7], one can construct the equilibrium measure corresponding to the Hahn polynomials. 
Recall c := n/N, and let 

a:=^-^Vl-c\ 6:=l + lyr^. (2.10) 

The density function is given by 

2 



arcsin( ^ ^ ^ ), x S (a, 6), 



TTC '2/ 



X — X 



(2.11) 

X € [0,a] U [b, 1]. 



1 

^ c 

Observe that the equilibrium measure of the Hahn polynomials corresponds to the saturated-band- 
saturated configuration defined in [1], and we also note that 

a + b = l, a-b = c^/4. (2.12) 
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2.2 The basic interpolation problem 



We begin with the basic interpolation problem for the Hahn polynomials. Note that the leading coefficient 

of PN,n{z) is 

iV"(n + a + /3 + l)„ 



(a + lU-N + l) 
Hence, the monic Hahn polynomials TTN,niz) are given by 



1 



k 



N,n 



PnAz)- 



(2.13) 



(2.14) 



Clearly, they satisfy the new orthogonality relation 

7V-1 



(2.15) 



fe=0 



where 7^^^ = k%^Jh*^^^. 

Following [1], we construct the interpolation problem for a 2 x 2 matrix- value function Y(z) with the 
properties: 

(Ya) Y{z) is analytic for z G C \ Xj\f; 

( Yfc) at each x^^k G ^N, the second column of Y has a simple pole where the residue is 

WN,k\ 



Res Y{z) = lim Y{z) , 

Z=XM,k Z-^XMh \ 



(2.16) 



( Yc) as z — )• 00, 



Y{z)=[l + 0[- 



z'' 




By the well-known theorem of Fokas, Its and Kitaev [5] (see also Baik et al. [1]), we have 
Theorem 2.1. The unique solution to the above interpolation problem is given by 



( 



Y{z) :-- 



T^N,n{xN,k)'^N,k \ 



k=0 ^ ~ ^N,k 



2 , , ^^^lN,n-l'^N,n-l{xN,k)wN,k 

rN,n~l^N,n~l{z) E 

A,-=0 



Z - XN,k 



3 Riemann-Hilbert Problem 



(2.17) 



In this section, we will make a sequence of transformations: Y ^ H ^ R, and transform the basic 
interpolation problem into an equivalent Riemann-Hilbert problem (RHP). 
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First, we want to remove the saturated regions. Following Baik et al. [1], we introduce the first 
transformation 



H{z) := Yiz) 



^ n - XNj)-' 

j=0 



V 









7V-1 

n - xN,j] 

j=0 



\ 



(3.1) 



It is readily verified that H{z) is a solution of the interpolation problem: 
( Ha) H{z) is analytic for 2 G C \ X^] 

{Hi,) at each XN,k £ ^A^) the second column of H{z) is analytic and the first column of H{z) has a simple 
pole where the residue is 

/ 



Res H{z) = lim H{z) 







1 



AT-l 



n {xN,k - XN,j) 



(3.2) 



( He) as z — oo, 



H(z) ={1 + 







Next, we shall employ an idea of Ou and Wong [10] to remove the poles in H{z), and transform this 
interpolation problem into an equivalent RHP; see also [15]. Let (5 > be a sufficiently small number, 
and we define 

/ 1 \ 



R{z) := H{z) 



for z G il.±, and 



-p^g±7V7rj2 cos{N7Tz) 
N-1 

Nttw{z) n iz~XNj)'^ 
j=0 



R{z) := H{z) 



(3.3) 



(3.4) 



for z ^ n±, where the domains are shown in Figure 1. Let S+ be the boundary of r2+ in the upper 
half-plane, and S_ be the mirror image of in the lower half-plane. For the contour S = (0, 1) U T,±, 
see also Figure 1. 

Lemma 3.1. For each xn k £ X^, the singularity of R{z) at xjy k is removable; that is, Res R{z) = 0. 



Proof. For z G ft±, it follows from (3.3) that 



Ruiz) = Hu{z) + Hi2iz) 



-p^g±7V7r»2 cos(A^7rz) 



Ruiz) = H^iz). 



(3.5) 
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l + i6 



1-iS 



Figure 1: The domains il.± and the contour S. 



Since Hi2{z) is analytic by (Hh), the residue of Ri2{z) at x^^k is zero. From (3.2), we also have 



Res Hii{z) = Hi2{xN,k) TT i^N,k - xn. 

'=XN,k WNk 



j=0 



On the other hand, it is readily seen that 



^^g_i_iv n (;;os(iV7rz) 



Res Hi2{z)- 



N-1 



Hl2{xN,k) {xN,k - XN,j) 



WN,k 



j=0 



(3.6) 



(3.7) 



Thus, applying (3.6) and (3.7) to (3.5) shows that the residue of Rii{z) at z = xjv,^ is zero. The entries in 
the second row of the matrix R{z) can be studied similarly. This completes the proof of the lemma. □ 

Note that this transformation makes R+{x) and R-{x) continuous on the interval (0,1). As a conse- 
quence, we have created several jump discontinuities on the contour S in the complex plane. 
It is easily verified that R{z) is a solution of the following RHP: 

(Ra) R{z) is analytic for z G C \ S; 

{Rb) the jump conditions on the contour S: for x G (0, 1), 

R+{x) = R-{x) 



1 0' 

r{x) 1 



(3.8) 



where 



for z £ E-t, 



r{x) 



icos'^iN'Kx) 



Ciw(x) llf=o^{x - XN,jY 



ci := 2N-Ki] 



R+{z) = R_{z) 



1 0^ 

?±{z) 1 



(3.9) 



(3.10) 



7 



where 



7±{z) = ' ^; (3.11) 



(Re) as z — oo, 



R(^) = {i + 0{- U ^ ^_„,„ L (3.12) 



4 The Auxiliary Functions 

To normalize the condition (Re) in the RHP for R{z), we need some auxihary fmictions. 
Definition 4.1. The g -function is the logarithmic potential defined by 

g{z):=[ log{z-s)fi{s)ds, zGC\(-oo,1], (4.1) 

Jo 

where iJ-{x) is the density function given in (2.11), and the so-called (p-function is defined by 

<P(z):=l/2-g{z), zeC\(-oo,l], (4.2) 

where I := 2 log \a — s\^{s)ds is called the Lagrange multiplier, a being the constant given in (2.10). 
Now, we define 

g{z):=giz)--^iz), (4.3) 
c 

where 

(p{z) = f log(z - s)ds. (4.4) 
Jo 

In (4.1) and (4.2), we view log(a; — s) as an analytic function of the variable s with a branch cut along 
(x,+oo). Since 

log_|_(a; — s) = log jx — s| lb iri, s G {x, +oo). 



we obtain 



PX pi 

9±{x) = / log{x — s)iJ,{s)ds -\- / {log \x — s\ ±TTi)^{s)ds 

Jo Jx 

= / \og\x — s\^{s)ds r^Tii / fi{s)ds. 

Jo Jx 



Similarly, we have 



/■X pi pi 

ip±{x) = / log \x — s\ds -\- / \og^{x — s)ds = / log \x — s\ds ± TTi{l — x). 

Jo Jx Jo 



8 



From (4.3), it follows that 

^ 1 1 

9±{^) = / log |x — s|(/i(s) )ds±'Ki / )ds, x G (0, 1) 

Jo c c 

Moreover, it can be easily verified that the ^-function satisfies the jump conditions 

1. 



a+{x) - g-{x) = 2m{l ), xG(-oo,0) 



(4.5) 



(4.6) 



and 



9+{x) - g-{x) 



X G (0, a), 



27ri(l - -), 

c 

2iri / (fJ,(s) )ds, x G (a,b), 

.L c 



(4.7) 



0, 



X G (6,1) 



On account of (4.1), (4.6) and (4.7), one readily sees that e"'^^^^ can be analytically extended to C \ [a, 6] 
and 



[1 + 0(i)] as Z — )■ OO. 



Hence, by adopting the convention that denotes the Pauli matrix 



c^3 



'1 ' 
-1 



we introduce the transformation 

T{z) := (cie"0-'^3/2^(^)g-ng(^)a3(g^gni^<73/2 



(4.8) 



(4.9) 



(4.10) 



to normalize the behavior of R{z) at z = oo, where ci is given in (3.9). 

As we have mentioned before, the density function /u(x) reaches its upper constraint near the endpoints 
of the interval of the orthogonality; see (2.11). Furthermore, since fJ.{x) is not differentiable at the point 
X = a, the function (j){x) is not analytic in the neighborhood of x = a, and we can not construct our 
global parametrix (such as Airy parametrix) by using 4){z) in the region which includes the critical point 
a. Hence, for our future analysis, a few more auxiliary functions are needed. 

Note by (4.7) that g{z) is analytic for z G (6,1). This together with (4.5) evokes us to consider a 
modified measure fJ-*{x) = fi{x) — i. Furthermore, we want to find a complex function z/(z) which satisfies 
the requirement 



t^±ix) = ±7ri(/i(x) 



1. 



(4.11) 



for x G (a, b). 
Let 



uiz) := - log[c/2 - y/{z-a){z-b)] - - log(^ - z^), 



(4.12) 
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where ^/{z^^a){z'^^V) is analytic in C \ [a, b] and behaves Hke z as z ^ +00, and where log( 
analytic function with branch cuts along (— oo,0] U [l,+oo). 
In view of the identities 



z — z'^) is an 



arccos x 



■ilog{x + i\/l — x^), — arccos x = —i log{x — i\J 1 — x^) (4.13) 
for X G [—1,1], we have 

2 c 

v±{x) = ^-iwcco^i^ — )■ xG(a,6). (4-14) 

c 1\Jx - 2;2 

On the other hand, the measure \i{x) — 'i/c can be rewritten as 

12 c 
— - = — — arccos( ^ ^ ^ i), x ^ {a,h). (4-15) 



c vrc 2V; 



, a; — X 

Prom (4.14) and (4.15), it clearly follows that 

u±{x) = ±-Ki{^(x) --), xe(a,b). (4.16) 

c 

Now, we are ready to introduce the auxiliary function 

i^{s)ds, z G C\ (-00,0] U [a,+oo), (4.17) 



where the path of integration from a to ^ lies entirely in the region ^ G C \ (—00, 0] U [a, +00), except for 
the initial point a. 
Similarly, we set 

(j)*(z):=l iy{s)ds, z G C \ (-00, 6] U [1, +00), (4.18) 
Jb 

where the path of integration from b to z lies entirely in the region z £ C \ (—00, b] U [1, +00), except for 
the initial point b. 

The functions (piz) and 4>*{z) defined above play an important role in our argument, and the following 
are some of their properties. 

Proposition 4.2. The mapping properties of <j){z) are given by: 

Re(/)±(x) < 0(0), Im0±(x) = =F-7rx, xG(-oo,0); (4.19) 

c 

0(x)<O, Im0(x) = O, xG[0,a); (4.20) 

^(a) = 0, ^±(6) = ±(1 - -)7ri; (4.21) 

c 

Re^±(x)=0, arg(^±(x) =^^, x e {a,b); (4.22) 

Re^±(x)<0, Im^±(x) = ±(1 - -)7r, xG(6,1]; (4.23) 

c 

Re^±(x) < Re0±(l), Im^±(x) = ±(1 - -x)7r, xG(1,+oo). (4.24) 

c 
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Proof. From the definition in (4.12), we note that ^{z) is analytic in C \ (—00, 0] U [a, b] U [1, +00). When 
X G (—00,0), it is readily seen that 

= z^i(x) =F -vri, (4.25) 
c 

where 

ui{x) = - log[c/2 + y/{a-x){b-x)] - - log{x'^ - x). (4.26) 
From (4.17) and (4.25), we have 

4'±{x) = / i'{s)ds + / i'±{s)ds 

J a Jo 

= ^(0)+/ ui{s)ds^-7rxi. (4.27) 
Jo c 

It is easily verified that z^i(s) > for s E (— oo,0), and hence (4.19) holds by (4.27). 

For X £ [0,a), i'(x) is analytic and it can be shown from (4.12) that i^(x) is real and positive. Thus, 
from (4.17) it follows that 

lm^(x)=0 and ^{x) < 0; (4.28) 

i.e., (4.20) holds. 

On the other hand, (^(a) = follows immediately from the definition in (4.17). Note that i'±{x) = 
ib7ri(/i(x) — i) for x G {o.,b); see (4.16). Hence, (4.21) holds by (4.17). If x G {cL,b), we also have from 
(4.16) 

4'±{x) = ±TTi[ {n{s)--)ds. (4.29) 

Ja C 

Note that /i(s) - ^ < for s G (a, 6) by (2.11), and hence (4.29) gives (4.22). 
For X G (6, 1], it follows from (4.17) that 

4'±{x) = / v±{s)ds + / i'{s)ds 

J a Jb 

= ^±(6)+ / iy{s)ds. (4.30) 



By a change of variable s = 1 — it can be verified that 



l-x 



v{s)ds= / -u{l-C)dC 
1 

-l-x 



(j){l-x). (4.31) 
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0(0) 



0(a) 





</>+(!) 



(a) The upper half z-plane under (j){z) (b) The lower half z-plane under (p{z) 

Figure 2: The image of the z-plane under the mapping (p{z). 



Applying (4.21) and (4.31) to (4.30) yields 



^±{x) = 0(1 -X) ± (1 - -)TTi. 

c 



(4.32) 



Thus, (4.23) holds by (4.20) and (4.32). 

For X £ {1, +oo), we have 



u±{x) = U2ix) =F -vri, 
c 



where 



'^2ix) = - log[y/ (x - a)(x - b) - c/2] log(. 



X — X 



(4.33) 
(4.34) 



From (4.23) and (4.33), it follows that 

0±(x) = 



1 PX 

iy±{s)ds + / iy±{s)ds 

'/'±(l)+/ iy2{s)dsT -{x -l)Tri 
Ji c 

Re0±(l)+ / i/2(s)ds ± (1 - -x)m. 
Ji c 



Note that 1^2(3) < for s G (1, +00). Thus, iy2{s)ds < and (4.24) holds. 



(4.35) 



□ 



A geometric interpretation of the results established above is given in Figure 2. Furthermore, if we 
let C+ = {z : Imz > 0} and C_ = {z : Imz < 0}, then for z £ C± we have from (4.16) and (4.17) 

(l)*(z)=f v{s)ds= I v{s)ds±-Ki I ( ^(s) - - ) = 0(z) ± 7ri(- - 1). (4.36) 

Jb Ja Jb \ Cj C 

The results in the following proposition can now be proved either directly from definition (4.18), or by 
using (4.36) and Proposition 4.2. 
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Proposition 4.3. The mapping properties of <j)*{z) are given by: 

Re (x) < (/>*(!), Im(^^(x) = ±^(1 -x)7r, xG(1,+oo); (4.37) 

^*{x)<0, lm^*{x)=0, xe{b,l]] (4.38) 

(j)*(b) = 0, (l)X(a) = ±(- - l)7ri; (4.39) 

c 

Re(/)^(x)=0, arg(/)i(x) = x £ {a,b); (4.40) 

Re(^±(x)<0, Im(/)i(x) = ±(^ - l)7r, xG[0,a); (4.41) 

Re(/>±(x) < Re(/.±(0), lm(l)*^{x) = ±{^—^ - -x)7t, xG(-oo,0). (4.42) 



The images of the upper (lower) half of z-plane under the mapping (j)*{z) is depicted in Figure 3. The 
next proposition shows the relationship between the function and the 0- function ((^* -function). 

Proposition 4.4. With (j){z) defined in (4-2), the following connection formulas between the (j)-function 
and the (j)- function ((f)* -function) hold: 

^{z) = 4){z) ± 7ri(l - -z) (4.43) 
c 

and 

(t)*iz) = (t)iz)±m-{l- z) (4.44) 
c 

for z G C-t . 




(a) The upper half z-plane under (t)*{z) (b) The lower half z-plane under (t)*{z) 

Figure 3: The image fo the z-plane under the mapping (t)*{z). 
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To establish (4.43) and (4.44), we need some properties of g'{z). In fact, we shall derive explicit 
formulas for g[z) and g'{z). From (2.11), we have 



- / log(2: — s)ds H — / log(z — s)ds + / log(z — s)fj.{s)ds. (4.45) 

C Jo C Jh J a 



Integration by parts gives 



log(z — s)fi{s)ds = /i(s) log(z — s)s — / s[log{z — s)fi{s)]'ds 



fi{s)log{z — s)s — / iJ,{s)ds — / slog{z — s)iJ,{syds + z / fi{s)ds. 

Ja J a Ja Z — S 

(4.46) 



The last integral in (4.46) is given by 

rb 



f 1 ^ f I 

/ ^{s)ds = —^{s)\og{z — s) + / ^{s)' \og{z — s)ds. 

Ja Z - S Ja 



On account of (2.11) and (2.12), a direct calculation gives 



and 



ix{x)' 



fi[s)ds = 1 

c 



1 - 2x 



2vr (x - x'^)'s/{x - a){b - x) ' 
Applying (4.49) to (4.47), we obtain 



1 



1 



z — s 



Similarly, one can show that 

rb 



fi{s)ds = — log(z - 6) + - log(z - a) - — 



X G (a, b). 



1 /"^ log(z - s) 



c c 

rb 



1 



\og{z - s) 



s \og{z — s)fi{syds 



27rA (l-s)V(s-a)(6-s) 
log(.-.) 



2vrA sy/{s-a){b-s) 
ds. 



ds 



\og{z - s) 



T^Ja ^{s-a){b-s) 27rA {I - s)^ {s - a){b - s) 
A combination of (4.46), (4.48), (4.50) and (4.51) gives 



log(z — s)fi{s)ds 



-(z — a) log(z — a) (z — b) log(z — b) 

c c 



z 



log(z -s) , 1 . 

.ds + —{z-l) 



\og{z - s) 



(4.47) 



(4.48) 



(4.49) 



(4.50) 



ds. (4.51) 



1 

+ - 



2TTJa s^{s-a){b-s) 2Tr' ' (1 - s) V(s - a)(6 - s) 
log(z - s) 

T^Ja ^J{s-a){b-s) 



zds 



ds. 



(4.52) 
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Let the integrals on the right-hand side of the equahty be denoted by Ii, I2 and I^, respectively. To 
evaluate Ii, we note that from [16] we have 

rb 



log s 



zds 



2tt 



log 



z + ^/^+(z-a)i/2(z- 6)1/2 



la {s-z)y^{s-a){b-s) (z-a)i/2(^_ 6)1/2 
Making the change of variable s = z — x, we obtain from (4.53) 

log(z - s) 



(4.53) 



a Sy^ {s — a){h — s) 



zds 



logx 



-a [z - x)y' {z - X — a){h - z + x) 

47r , z + J{z-a){z-h) + c/2 

log 

c [sjz — a + \/ z — b)z 



d{z — x) 



(4.54) 



and 



log(z - s) 



a (l-s)V(s-a)(6-s) 



zds 



47r , Z-1 + J(z - a)(z - b) - c/2 

log z^=^= . 

c {y/z - a + \/z - b){z - 1) 



Also, note that 



/ 

J a 



log S 



see [16, Lemma 1]. Hence, 

^3 



V^(s - a){b - s) 
' log{z - s' 



ds = 27rlog 



y/E+Vb 



a ^/{s^^a)ib^^ 
A straightforward calculation shows that 



ds = 27r log 



\J z — b + \J z — a 



(4.55) 



(4.56) 



(4.57) 



g[z) = -1 - 21og2 + -(z - l)log(z- 1) - -zlogz+ (2 - -)log \[z - a)^/2 + (^ - hfl'^ 
c c c V 

+ -zlog U + (z - a)i/2(z - 6)1/2 ^ ^/2l + -(1 - z) log L - 1 + (z - afl'^^z - hfl'^ - c/2 
c L J c L 



(4.58) 



Moreover, by using (4.50), (4.54) and (4.55), we have 

g{z) = I ^^^^ ds — - log(z — a) + - log 2; — - log(2; — 1) + - log(z — 6) 
J a z — s c c c c 



2, 

-log 
c 



z + .J{z-a){z-b) + c/2 



{z-l) + ^/{z-a){z-b) -c/2 c 



I, fz-l 
+ -loi 



(4.59) 
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Proof of Proposition 4-4 ■ Note that when x S (a, 6), i'±{x) = ±7ri{ij,{x) — i); see (4.16). Hence, it follows 
that 

u+{x) -u^{x) = 2m{fi{x) --). (4.60) 

c 

On the other hand, if a; G (— oo, 0) U (1, +oo), we have with the aid of (4.25) and (4.33) 

2 

i^+(x) — i^-(x) = ni. (4-61) 

c 

To establish (4.43), we need one more function (j){z) defined by 

^{z) := -g'{z)=f^TTi, z e C±. (4.62) 
From (4.1), one can show that for x G (0, 1) 

g'+{x)-g'_{x) = -27rin{x). (4.63) 

From (4.63), it is easily verified that (p{z) is analytic for z G C \ (— oo, 0] U [a, b] U [1, +oo) and satisfies the 
jump conditions 

4>+(x) - 4>^(x) = 27ri(n(x) --), xe(a,b) (4.64) 

c 

and 

2 

(l)+{x) - (j)_{x) = --TTi, x E (-00,0) U (l,+oo). (4.65) 
Comparing (4.64) and (4.65) with (4.60) and (4.61), we find 

i^+{x) — z/_(x) = (j)+{x) — (j)-{x) (4.66) 

or, equivalently, 

i^+{x) — 4>+{x) = v^{x) — (j)-{x) (4-67) 

for X £ (—00,0) U {a,b) U (l,+oo). On account of this, z/(z) — (j){z) can be analytically extended to the 
whole complex plane. Note that i^+(x) = — j^_(x), and that from (4.59) it can be shown 

//N //N 2, / X + xc \ 2, /I — x\ , 

g;(x) + g^_(x) = -log(^ ^^_^^^^^^J +-log [—^) =0 x€{a,b) (4.68) 

and so 

4>^{x) = -^_{x), xe{a,b). (4.69) 
Hence, from (4.67), (4.16) and (4.69) we have 

u{x) - ^{x) = l[{i^- ^)+{x) + {iy- 4>)-ix)] 

= 0, x£{a,b). (4.70) 
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Using analytic continuation, it follows from (4.62) 

u{z) = -g'{z) =F -vrz, z e C±. (4.71) 
c 

Hence, it follows from (4.17) that 

4>{z) = j g'{s)^-^-Ki^ds, zgC±. (4.72) 
Now, let z G C+; on account of (4.1) and the definition of Z, we have 

(y - a{s) - ^-Ki^ds 

= -9{z) + g+{a) - -7Ti{z - a) 
c 

= 1/2- g(z)+m(l--z). (4.73) 
c 

In a similar manner, one can show that for z € C_ 

4){z) = l/2-g{z)-TTi{l-^z). (4.74) 

Formula (4.43) now follows from (4.73), (4.74) and (4.2). Formula (4.44) is obtained by coupling (4.43) 

and (4.36). □ 

We conclude this section with the mapping properties of (j){z) in (4.2). By using Proposition 4.2 and 
(4.43), one can prove the following proposition; see also Figure 4. (To establish the second equality in 
(4.78), we make use of (4.44) instead of (4.43).) 

Proposition 4.5. The mapping properties of <j){z) are given by: 

Re(j)±{x) <Re(j)±{0), Im (j)±{x) = ^ir , x £ {-oo,0); (4.75) 

Re(l)±(x)<0, lm(f)±(x) = ^(l--x)7T, xe[0,a); (4.76) 

c 

Ma) = T{1 - -a)TTi, Mb) = T-{l-b)7ri; (4.77) 
c c 

Re0±(x) = O, arg(/.±(x) = T^, x £ {a,b); (4.78) 

Re4>±(x)<0, Im(/>±(x) = =F-(1 -2;)7r, x e (b,l); (4.79) 

c 

(pix) < (pil) < 0, x€ (1, +oo). (4.80) 
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(a) The upper half z-plane under (t){z) (b) The lower half z-plane under (^{z) 

Figure 4: The image of the z-plane under the mapping (/>(z). 

5 Normalization and Contour Deformation 

Recah from (4.10) that 

T{z) = (cie"')-"^/2^(z)e-"^(")"^(cie"')"^/2^ (5.1) 

where / is the constant given in Definition 4.1 and g{z) is shown in (4.3). By using (3.8), (3.10), (4.8) and 
(3.12), a straightforward calculation shows that T(z) is the unique solution of the following RHP: 

(Tq) T{z) is analytic for 2; G C \ S; 

(Tb) for z £T.±, 



Mz) = T.{z) 



and for z € (0, 1), 



where 



Mz) = T^{x) 



1 0^ 

cif±(z)e"('"23W) 1^ 

Jll{x) Jl2{x)\ 



J2l(x) J22{x) 



(5.2) 



(5.3) 



Jii(x) = e-"(^+W-9-W), 
J2i(x) = cir(x)e"('-^-(^)-^+(^)), 



and 



Ji2(a;) = 0; 

(Tc) T{z) behaves like the identity matrix at infinity: 

T{z) = I + 0{l/z) 



as z — )• 00. 
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For the sake of simplicity, we introduce some auxiliary functions. Let 

„1 N 



E{z) := exp ^-iV j \og{z — s)ds^ 11 (-^ ~ ^Nj) 



N-l 



-Nip{z) 

j=0 

for z G C \ [0, 1] and 



E(z) := ^ Z^, , (5.5) 
^ ' 2cos{N7Tz) ^ ' 

for z € C±, where (p{z) is defined in (4.4). In view of the equation preceding (4.5), it is easily verified 
that 

E+{x)/E_{x) = e^^(^- W-v+W) = e2^™, x e (0, 1). (5.6) 

On account of this, it can be shown that E^{x)E^{x)"^ = 1 for x G (0, 1). Hence, E{z) can be analytically 
continued to the interval (0, 1). Moreover, 

4cos2(iV7ra;) 

As — 7- oo, by applying Stirling's formula to (5.4) and evaluating the integral in (4.4) explicitly, one can 
show that E[z) ~ 1 uniformly for z bounded away from the interval [0, 1]. Moreover, from (5.5), we have 
for z G C± 

E{z) /E{z) = 1 + e^'^^^'' ~ 1 (5.8) 

uniformly for z bounded away from the real line. 

To compute the jump conditions of T{z) in (T^) above, we first state the following lemma. 

Lemma 5.1. For < x < 1, we have 

r(a;)e"('-3+(^)-^'-(^» = ^ , (5.9) 

ciw{x)E{xY 

and 

_ / 1, rE G (0, a) U (6, 1) , 

" \ e2"<^-(-) = e-2-^;(-), X G (a, 6). ^^'^"^ 

For z G S±, we have 

p2n<j>{z) 

f±(z)e"('-23(-» = . (5.11) 

^ciw{z)E{z)E{z) 
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Proof. For < x < 1, it follows from (3.9) and (4.3) that 

4cos^(iV7rx) 



r{x)e' 



n[l-g+{x)-g-{x)) 



^n{l-g+-g-) 



ciwe-^^^++v-) W]~q {x - XN, 
4cos^(iV7ra;) ^„(^*^^*^ 



(5.12) 



The second equality is obtained by using (4.2), (4.44) and (5.4). Substituting (5.7) into (5.12) gives 
(5.9). For X G (0, a) U {b,l), the formula e^^^+^^^~^-^^'^^ = 1 follows immediately from (4.7). Note that 
i'±{x) = ±iTi{fi{x) — -) for X G {a,b); hence, by (4.7) and (4.18) we obtain 



^nig+{x)-g-(x)) _ ^2nJ^ u-{s)ds _ ^2n<j>t{x) _ ^-2n4>\{x) ^ ^ ^ (a,6). 

For z £ T,±, (5.11) follows directly from (4.2), (5.4) and (5.5). 



(5.13) 



□ 



With the properties of (p{z) and 4'*{z) established in Section 4 and Lemma 5.1, the jump matrices for 
T in condition (T;,) can be written as follows: 



r+(x) = r_(x) 



1 







^/{wE^) 1 



for < X < a. 



for a < X < b, 



for b < X < 1, and 



r+(x) = r„(x) 



r+(x) = r_(x) 



T+iz) = T_iz) 



l/{wE^) e 



-2ntl>*, 



1 ^ 0' 



1 



0' 



^e^'"l'/{wEE) 1 



(5.14) 



(5.15) 



(5.16) 



(5.17) 



for z G Tj±. 

Note that for the jump matrix on (a, 6), we have the following factorization: 



\l/{wE'^ 




(e wE/e^''^- 
\0 l/E 




(l/E wE/e^'"^*+ 



(5.18) 



V 



E 



The first and third matrices on the right-hand side of (5.18) have, respectively, the analytic continuations 



E wE/e^^'f'* 
l/E 



and 



l/E wE/e^'^'^* 
E 



on both sides of (a, 5). Let S-t = |J-^-|^S!|_, and let Sy 
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l^i^i U Sj. U denote the contour shown in Figure 5. Define 



Viz) :=T{z)E{zr, 
V{z) :=T{z)E{zr\ 



V{z) := T{z) 
V{z) := T{z) 



,0 l/E 

'e w^/e2"<^*' 
l/E 



z G 



(5.19) 
(5.20) 

(5.21) 
(5.22) 



where regions 0± = U^^^ fi*^ and n are also depicted in Figure 5. It can be verified that the jump 
matrices of V{z) are given by 



Jvix) 
Jv{x) 
Jv{x) 



-w\ 
) ' 

' 1 0^ 



z G (a, 6); 
z G (0,a); 
zG (6,1); 



(5.23) 
(5.24) 
(5.25) 



0. 



















\ E^ 






a \ 




h 










1 E^_ 





S3 



EL 



EL 



fin 



Figure 5: Contour Sy and regions ^j-i ) ^oo* 
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and 



Jv(x) 



E/E J) \ 



(5.26) 



Jv{z) 



I 1 ±wEl{e^'"i'* E)" 



w e 



E/E 



(5.27) 



(5.28) 



Jv{z) 



1 =Fw/e 




(5.29) 



In proving (5.23), (5.24) and (5.25), we note that E is analytic in < x < 1. For (5.23), we also make 



use of (5.18). In establishing (5.27), we have used (4.44) and (5.8) 
T_(z)^^T+(z) = /. Thus, V{z) is a solution of the following RHP: 



For (5.28) and (5.29), we note that 



(K) V{z) is analytic for z E C \ Sy U [0, 1]; 
(T4) for z G Sy, 

(Vc) as z — oo. 



V+iz) = V.iz)Jviz); 



V{z) = I + 0{l/z). 



For condition (Vc), we note from (5.4) that E{z) 1 as z — )• oo. From the properties of (l)*{z) and <l){z) 
(see also Figure 2 and Figure 4), the regions and can be chosen sufficiently small so that 



Re<p*{z)>0 for z G 
Re(/)(z) < for z G il.±. 



(5.30) 



Moreover, from (5.8) it follows that E/E ~ 1 as n — t- oo for z bounded away from the real line. These 
together with (4.20) and (4.38) imply that the jump matrix Jyiz) tends exponentially to the identity 
matrix as n — )■ oo for z bounded away from [a, 6] U {0} U {!}. On the other hand, the weight function 
w{z) in (1.7) can be rewritten as 



w{z) 



a\l3\N^°^^l^ p{Nz - l/2;a,p,N - 1) 
TjNz + a + l/2)r(jV(l -z)+/3 + 1/2) 
N'^+PT{Nz + 1/2)T{N{1 - z) + 1/2) 



(5.31) 
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A direct calculation shows that 

w{z) ~ - zf as iV ^ oo. (5.32) 

Therefore, it is natural to suggest that for large n, the solution of the RHP for V{z) may behave asymp- 
totically like the solution of the following RHP for N{z): 



{Na) N(z) is analytic for z G C \ [a, 6]; 
(Nb) for X G {a,b), 



A.,W = I ° -'f 1, ,5.33) 



where 



(Nc) as z — )• oo, 



h{x) = x'^{l-xf; (5.34) 



N{z)=I + 0{l/z). 
To solve the above RHP, we first consider the function 



/3 



. z + c/2+{z-ay/\z-b)'/' \ l-z + c/2-{z-ay/Hz-by/' 

V z{^ + Vb) J \ {i-z){^ + Vb) J ' 

where (z — a)^^'^{z — 6)^/^ is analytic in C \ [a,b] and behaves like z as z — )■ oo. In view of (2.12), it is 
easily shown that the function M{z) is indeed analytic in C \ [a, b] and 

M+(x)M_(x) = x""(l - x)-/^ for x G (a, b). (5.36) 

Its limit at infinity is given by 

Moo = lim M{z) = (2/{^+Vb)Y^^ . (5.37) 



2— >00 



Therefore, if we set 

N{z) := M^N{z)M{zy^, (5.38) 

it follows from (5.33), (5.36) and (5.37) that N{z) is a solution of the RHP: 
(Na) N{z) is analytic for z £ C\[a,b]; 
(Nb) for X G (a, 6), 

/V^j^r.Tl = /V_f.Tl I 

1 



N+{x) = N^ix)r (5.39) 
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{Nf.) as z — )• oo, 



N{z)=I + 0{l/z). 



This problem can be solved explicitly, and its solution is given by 



N(z) = -( «(^)+«(^)"' i{aiz) - a{z)"^) 
^ ' 2 U(a(z)"i -a(z)) a(z) +a(z)-i 



where 



a{z) 



1/4 



(z-a)V4 



(5.40) 



(5.41) 



with a branch cut along [a, 6] and a(z) — t- 1 as z — t- oo. 
Hence, 



iV(z) 



a(z) + a{z) ^ i(a{z) — a(z) ^) 



2"^^ \i{a{z)-^ - a{z)) a{z) + a{z)-^ 
and for z G C \ {0, 1, a, b}, we expect that 

v{z)r^N{zy, 



M{z] 



0-3 



(5.42) 



(5.43) 



see the comment following equation (5.32). However, since E/E 9^ 1 for z near and 1 and since 
g-2n</>*(2) ^ g £qj, ^ near a and 6, the jump matrix Jv{z) 7^ I near these critical points; see (5.26)-(5.29). 
That is, V{z) and N{z) are not uniformly close to each other near these points, and special attention 
must be paid to the neighborhoods of these points. 

Let U (0, e) := {z £ C : \z\ < e} be a neighborhood of the origin, where e is a small positive number. In 
view of (5.5), (5.26), (5.30) and (5.43), we first wish to find a matrix-valued function Vo{z) that satisfies 
the following RHP: 

(^,a) yo{z) is analytic for z £ U{0,e) \ {—i£,ie); 



(VoUiz) = {Vo)-{z) { 



1 + e^^'^*^ 







(1+ 



for z G (0, +ie), 
for z G {—is, 0), 



(5.44) 



/ 1 + e-'^^'^'^ 
1^ (1 + e-2^'^*^)-iy 

where the functions {Vo)-^-{z) and (Vo)_(z) denote the boundary values of Vo{z) taken from the left 
and right of the imaginary axis, respectively; 



{Vo,c) fovz£dU{0,e) 



Voiz) ~ N{z) 



(5.45) 



as n — )• 00. 
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Next, we consider as in [15] the following D-function 

^ e^^T{Nz + 1/2) 

By applying Stirling's formula to (5.46), we have 

D{z) = 1 + 0(l/n) 
for I arg^l < vr. We also introduce the function 



z G C\ (-00,0]. 



as n — oo. 



D*{z) 



(l + e-2^'^^^)£)(z), zGC_, 



which is actually analytic in (— oo,0). From Euler's reflection formula, it follows that 

'/2^e^'{-Nz)~^' 



D*{z) 



T{-Nz + 1/2) 



zeC\ [0,+oo). 



Similarly, for | arg(— z)| < vr, we have 



D*iz) = l + 0{l/n) 



as n — )• oo. 

Set the following matrix-valued function 



Qoiz) = < 



D(z) 



, Diz)-') 

'd*{z) 
D*{z)~^ 



D{z)''\ Rez>0, 



L>*(z)<^3, Rez<0. 



It follows from (5.47), (5.48) and (5.50) that Qo is a solution of the RHP: 

iQo,a) Qoiz) is analytic for z G C \ (— ioo,«oo); 

{Qo,b) 

(1 + e2JV^*.)-i I fo'^ ^ e (0' 
'l + e-^^'^*^ 



{Qo)+iz) = {Qo)~iz) { 







[1 + e 



-2NiTiz\ 



-y for z G (— ioo, 0), 



(5.46) 



(5.47) 



(5.48) 



(5.49) 



(5.50) 



(5.51) 



(5.52) 



where the functions ((5o)+(-z) and ((5o)-(^) denote the boundary values of Qq{z) taken from the 
left and right of the imaginary axis, respectively; 



25 



{Qo,c) for z bounded away from the origin, 



Qo{z) ~ / 



(5.53) 



as n — )• oo. 



A comparison of conditions (Qo.fe) and (Qcc) with (5.44) and (5.45) shows that a solution of the RHP for 
Vo is given by 



Vo{z) = Niz)Qoiz) 



N{z)D{zy^, Re^>0, 
N{z)D*{zy-\ Rez<0. 



(5.54) 



In a similar manner, for z G U{l,e) = {z £ C : \z — 1\ < e}, we want to find a matrix-valued function 
Vi{z) that satisfies the RHP: 

{Vi,a) is analytic for z G U{1, e)\{l — ie, 1 + ie); 



(V^Uiz) = iV.Uz) 



1 + e^^'^*^ 

Q ^^ \ „2NTTiz\ — l 



l + e 



-2Nniz 











(l + e- 



-2NTTiz\ 



for z e {1,1 + ie). 



for 2: G (1 — ie, 1) 



(5.55) 



where the functions {Vi)^{z) and (Vi)_(z) denote the boundary values of Vi{z) taken from the right 
and left of the line Re^; = 1, respectively; also see Figure 5; 



(yi,e) lorz£dUil,e) 



Vi{z) ~ N{z) 



(5.56) 



as n — 7- 00. 

As in the case of Vq, it can be verified that a solution of the RHP for Vi is given by 



Viiz) 



N{z)D*{l - z)''\ Rez>l, 
N{z)D{l - z)"-^, Rez<l. 



(5.57) 



Since Vo{z) satisfies the same jump conditions as V{z) on the contour dU{0,e) U (0,+ie) U (— ie,0), 
V{z) ~ Vo{z) for z G U{0,e). Similarly, V{z) ~ Vi{z) for z G U{l,e). Hence, as n — )■ 00, we have from 
(5.43) 



N{z), z G (f^oo U U n^) \ {U{0, e) U U{1, e)) , 
V{z)r^{ Voiz), zeU{0,e), 
Viiz), zGU{l,e). 



(5.58) 
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By (5.19)-(5.20), it is clear that 



Tiz) 



v{z)E{z)~''-^, zGn^unl. 



(5.59) 



Let xo be an arbitrary fixed point in the interval {a,b). Prom (5.54), (5.58) and (5.59), it follows that 
for Kez < xq, 



T{z) ~ < 



z G Ooo\C/(0,e) 
zGJl|\C/(0,e), 



N{z)E{zy''--i, 
N{z)E{z)-''-^, 

N{z)D*{zY-^E{z)-"\ z G J7(0,e) nOoo, 

^ N{z)D{zY-^E{zy-\ z eU{^,e)r\Q%, 



(5.60) 



and by (5.47) and (5.50) we obtain 



T{z) ~ < 



' N{z)D*{zY^E{z)-''^, zen^\U{0,e), 

N{z)D{zY^E{z)-''\ z G \C/(0,e), 

N{z)D*{zY''E{zy"--^, z G ;7(0,e) nJloo, 

^ N{z)D{zY^E{z)-''^, z G ;7(0,e) nO|, 



(5.61) 



for Rez < xq. From the definitions of E{z) and D*{z) in (5.5) and (5.48), respectively, it is easily shown 
that 



D*{z)E{z) = D{z)E{z). 

Hence, it follows that for Rez < xq, 

T{z) ~ N{z)D*{zY^E{zy\ 
Similarly, for Rez > xq, we have by using (5.57), (5.58) and (5.59) 

T{z) ~ N{z)D*{l - zY-''E{zy\ z G J^oo U Jl^, 



z G Ooo U 



(5.62) 



(5.63) 



(5.64) 



where we have also made use of the identity D*{1 — z)E{z) = D{1 — z)E{z) obtained from (5.5) and 
(5.48). Define 



D{z) 



D*{z), Re z < Xq, 
D*{l-z), Re2;>xo. 



(5.65) 



On account of (4.3), (4.10) and (5.4), one can work backwards to get 

R{z) ~ (cie"')^='/2^(z)5(z)'^«e"^(^)'"«E(z)-'"^(cie"')-'"^/2 



(cie"')'^^/2^(z)e-"<^(^)° 



Diz) 



l/2T-r7V-l/ N 

Cl nj=0 [Z-XN,j) 



0-3 



(5.66) 



for z G r^oo U U O^, where / is given in Definition 4.1. 
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Figure 6: The regions I, II, III. 



6 Construction of the Parametrices 



In this section, we wih construct an approximation R{z) to the solution of the RHP for R{z) stated in 
Section 3 for large n. Given j;i G (0,a), define the rectangle -fC = {z G C : xi < Re^; < 1 — xi, |Imz| < 5}, 
and the line Fq = G C : Re 2: = xq, |Im z| < 5}, where xq is given in Section 5. We divide the complex 
plane into three regions I, II, III by using Fq and K; see Figure 6. 
From (5.66), it is natural to set 

~ 10-3 



D{z) 



Z — X]\fj J 



(6.1) 



for z £ I. 

For z inside the rectangle K, we will construct the parametrices by using Airy functions as in [3, 15]. 
For z G II, it follows from (5.48), (5.65) and (5.66) 

2cos{NTrz)D{z) "i '"^ 



R{z) ~ (cie"')'^^/2^(z)e~"'^(^)'^3e^^^^^'^^ 
and from (4.43) and (5.47) we obtain 



R{z) ~ (-I)"(cie"')'^3/2^(^)g-n0(^)a3 



.j=0 [Z - XN,j) 



2cos(7V7rz) 



z G C±, 



Cil/2 Y\f=0^{z-XN,j) 



0-3 



Using (4.12) and (4.17), one can derive an expansion for (p{z) for z G U{a,e) := {z G 
Indeed, we have 



3c2 



3/2 



(6.2) 

(6.3) 
Iz — a| < e}. 

(6.4) 



for z G U{a, e). By using (4.29) and the proposition of (j) (see Figure 2), it can be verified that the function 
defined by 

2/3 



(6.5) 
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satisfies 

f.{x)-'U{x) = l, xe{a,b). 

On account of this and (6.4), we conclude that this function is analytic in C \ (— oo,0] U [6, +00). In 
particular, one can see that for z G II n C-|_ 

-7r<arg7(z) <0, (6.6) 

and we obtain from (4.22) 

U{xf'^J.{x)-^'^ = e-^'l\ X G (a, 6). (6.7) 
Note that N{z) = M"''^ N{z)M{zy^ by (5.40) and (5.42), and N{z) has the factorization 



(6.8) 



Since a{z) = (z — bY^^ /{z — a)^/'^, ^{z) has fourth-root singularities at z = a and z = b. In order to 
remove the singularity at z = a, we introduce the function 

En{z) = N{z)h{zrl^ [n^l^J{z)r"^/\ (6.9) 

Note that by using (2.12), we have from (5.35) 

limM(z) = a-°/25-^/2. 

Also, from (5.34) we have 

lim h{z) = a"b^. 

Thus, 

M{zY^h{zY^''^ ^ I asz^a. 
Coupling (6.9) and (5.38), we obtain 

En{z) ~ M~^-^ ^. (ni/6G(z)/(z)^/4)~"' as z ^ a. (6.10) 

Since f{z) has a simple zero at 2; = a and a{z) has a fourth-root singularity at a, the last formula shows 
that En{z) has a removable singularity at 2; = a. 

Furthermore, one can show that En{z) has no jump across the interval {a,b). Indeed, by (5.33) 

N^{x)h{xr/' = N^{x) (^^^^°_, -^^^^^ h{xr/' 

= N^{x)h{xr/^ r xe{a,b). (6.11) 
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From (6.7), (6.9) and (6.11), it follows that for x G {a,b) 

1 -A /O -l\ ( 1 i 



1 



?l2/37_(^)]-3/4 



\-i 1 

-i o\ 

i / 



Ml / I i 1 / 

n2/37+(x)]-3/4 ^ ^. 



[n2/V4x)]-'^3/4 



(6.12) 



Therefore, En{z) is analytic across (a, 6). Also, note that from (5.42), (5.34) and (6.5), it is evident that 
N{z), h{z) and f{z) are analytic for z G (0,a). Thus, the matrix-valued function En{z) is analytic in 
C\ (-00,0] U [b, +oo). 

By virtue of (6.9), we can rewrite (6.3) as 



R{z) ~ (-l)"(cie"')"3/2^„,(^)^^(^)-i^(^)g-n0(2),.3 
(-1) 



2cos(iV7rz) 



CiV2 l\f^Q^{z - XN,j] 



0-3 



{c,e-Y'^'Er,{z)[n'/^fiz)r/' " ~\ "\ 



2cos(7V7rz) 



(ci/i(2;))2 nj=o - 



0-3 



To find an approximation to R{z), we first look for a matrix which is asymptotic to 

-n(p{z) 



n 



,3/4 e -'-v^^ -le-^y- 



Set 



e := n2/3/(z) = (- 
From the asymptotic expansions of the Airy function [9, p. 198], we have 



3 ^ \2/3 
--n(t){z) 



(6.13) 



(6.14) 



(6.15) 



Ai(0 



2^ 2^ 



Ai'(e) ~ -|^e-|«^^^ 



20F 



[n2/37(z)]i/4en^{.)^ 



2^ 



(6.16) 



2^ 



2^ 

i7r/6 



2/37(z)]"i/4e-"^(^), 



2^^^ 2^ 
where w = e^'^*/^ . It is immediate that for z G II n C+ , 



n2/3/(z)]i/4e-"^W, 



[„2/37(,)].3/4 [ ~ 2^ ^-^^'Ai'(^0 ^Ai'(e)^ 



-ze 



-n4>{z) pn4>{z) 



-LoAiicoo Ai(e) y 



(6.17) 
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Similarly, for z G II n C_, we also have < arg f{z) < vr and 



For z S III, there is a result corresponding to (6.3) with <j){z) replaced by (j)*{z). Indeed, it follows 
from (5.65), (5.66) and (4.44) that 



R{z) ~ (-l)^(cie"')'^3/2^(z)e-"<^*(^)^^ 



2cos(iV7rz) 



Let 



2/3 



which is analytic in C \ (— oo, a] U [1, +oo). By Proposition 4.3, we see that for x G (a, b) 
Since N(z) has a fourth-root singularity at z = b, we define 



(6.19) 



(6.20) 



(6.21) 



K(z) = N{z)hizr/' {\ ; 1 [n2/3r(^)]-3/4 



i 1 



(6.22) 



As in the case of En{z), it is readily seen that 5 is a removable singularity of E^{z). By (6.11) and (6.21), 
one also has 



(6.23) 



for X € {a,b). Therefore, E^{z) is analytic in C \ (— oo,a] U [l,+oo). On account of this and (6.19), we 
have for z G III 



R{z) ~ (-l)^(cie"')"--'/2£;*(^)^*(2)-l^(^)g-nrW-3 



2cos(A^7rz) 



0"3 



2 cos(iV7rz) 



(ciM^))i/2n;=o(^-^^,.) 



0-3 



Hence, as before, it can be shown that the matrix 

[n2/3/*(z)]-3/4 ( 



g-n</)*(2:) j^^n4i*(z) 
j^^-ncj)* {z) ^n(j>'{z) 



(6.24) 



(6.25) 
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is the leading term of the asymptotic expansion of the matrices 



'^^Ua;2Ai(a;2e) Ai(r) 



for z G Illn C_L, and 



2^ 



-lOJ 



2Ai'(wr) -iAi'(r) 



u;Ai(a;r) Ai(e) / 
for z e Illn C_, where C = ri?^'^f*{z) = (-|n(/)*(z))^/^. Define the matrix function 



P{z) :-- 



—UJ 



2Ai'(a^0 ^Ai'(J)\ 
Ai(a;e) Ai(0 ) 



.;2Ai(^2^) Ai(0 J ' 



-uj'^k:\{u?i*) Ai(e) 

a;Ai(wr) Ai(r) 



z G iinc+, 



z e iinc_, 



z G IIInC+, 



, zGiiinC- 



An appeal to the formulas 



shows P{z^ satisfies 



Ai(2;) + ujK\{ijjz) + bP'Kiiijp'z) = 0, 
Ai'(z) + Lo'^AS'iujz) + wAi'(a;2z) = 



P+(x) = P„(x) 



1 0' 



X G 



(6.26) 



(6.27) 



(6.28) 



(6.29) 



(6.30) 



Let En{z) and E^{z) be defined as in (6.9) and (6.22), respectively. In view of the heuristic argument 
leading to (6.1), (6.13) and (6.24), it is reasonable to suggest that R{z) is asymptotically approximated 
by 

' (cie"')'^"'^A^(^)e"'"^(^)''^Mi(z), z£l, 



R{z) := { {-lY^{cie^^Y^I^En{z)P{z)M2{z), z G II, 
^ {-l)^^{cie''Y'''^K{z)P{z)M2{z), zGlII, 



(6.31) 



where 



Mi(z) := 



D{z) 



Cl^/^nj=0 {z-XNj) 



0"3 



(6.32) 
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for z S I, and 



M2(z) : = 



2cos(iV7rz) 



0-3 



(6.33) 



for z Gil and III. From (6.31), it can be shown that R{z) has the same large z behavior as R{z) given in 
(Re)- To see this, we only consider the case when z £l and Kez > xq. The discussion for the case z & l 
and Rez < xq is very similar. Note that when z G I and Rez > xq, 



R{z) = (cie"')'"^/2^(z)e-"'^(^)'^3 



D{z) 



Cil/2 W^^Q{z-XN,j) 



(6.34) 



Prom (5.65) and (5.50), it is easily seen that D[z) = D*{1 — z) ~ 1 for large z. Furthermore, from (4.2) 
and rijL'o^l-^ ~ ^Nj) ^ z^ , we have 



R{z) ~ (cie"')^='/2^(z)(cie'^')-'"^/2gn3(^)a3^-7Va3_ 
This together with (4.1) and the asymptotic behavior of N{z) in {N^.) gives 

R{z) ~ (cie"')'^=*/2^(z)(cie"')"'''/^e"^'^-'^'''^"^''' 
= (cie'^')"^/2[j ^ 0(l/z)](cie"')-"^/2^("-^)"3 



(6.35) 



as z — 00. 



(6.36) 



Let r := IJ^^o ^« (0) 1); see Figure 7. Since R{z) has the same large z-behavior as R{z), it can be 
readily verified that the matrix-valued function 



5(z) := (cie"')-"^/2^(z)/?(z)-i(cie"')<^^/2 

is a solution of the RHP: 

(Sa) S{z) is analytic for z G C \ P; 

(Sb) for z G P, 



(6.37) 



where 



(5c) for z G C \ P, 



S+{z) = S.iz)Jsiz), 
Js{z) := (cie"0-'^^/2i?_(z)jR(z)i?+(z)-i(cie"')-^/2; 



(6.38) 



S{z) — )• / as z — >■ 00, 



(6.39) 
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Figure 7: Contour F. 



where the contour associated with the matrix Js{z) is F = S U Fq U dK; see Figure 1 and Figure 6. 

To solve this problem, we first derive the asymptotic behavior of the jump matrix Js{z) as n — ?• 00. 
For convenience, we consider only the cases when z lies on the contour Fq H C+ and F3 n C+, and z in 
the interval (0, 1), since the other cases can be discussed in a similar manner. 

For z G Fo n C+, we have from (6.31), (6.28) and (6.29) 



{-irV^{cie^Y'^^En{z) 



M2{z). (6.40) 



Similarly, 



R.(z) 



f^(r ."'^3/2^*/.^ f ^^Ai'(o;2e) zu;2Ai'(a;r)\ / 1 A , . . . 
-1) V^(cie ) /^n(^)l_^2Ai(^2^*) -^Ai(^e)JU l]^^'^'^- 



(6.41) 



Since -R(2;) is analytic in C\S, Jr{z) = 1 for z G FonC+; see the RHP for R{z) in Section 3. Thus, from 
(6.38) it follows that 

( ,.N-n^*,J ia;2Ai'(a;r)\ /-ia;'Ai'(a;0 -^a;Ai'(a;2^)\ ~ ^ 

\ -w Ai(a;^4 ) — a;Ai(a;4 ) / \ — a;Ai(a;0 — w Ai(a; 4) / 



On account of the well-known formula [9, p. 194, (9.2.9)] 



wAi(w^nAi'(a;0 - w^Ai(wOAi'(w^n = , 

27rz 



(6.43) 



we obtain 



^ ^ iuA\'iu^C) iu^A^uC)] -^'Aiiu^O i^Ai'(c^2^) \ ~ 1 

\ — 1^ Ai((^^^ ) -(^Ai((^4 ) / \ LoA\{ujk,) —luj'^ Ax [ojQ j 
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Now, we set out to derive the asymptotic behavior of Jsiz) for z S ronC4.. Since by (4.22) and (6.5) 
arg/(xo) = — vr, arg^ = arg/(2;) € [— vr, — vr/3) for z £ Tq nC+. We recah the asymptotic behavior of the 
Airy function and its derivative in [9, p. 198], and have as in (6.16) 

-w2Ai(a;2^) = |-i^e"^(")(l + 0(l/n)), iujAi'{oj^^) = ^e"^W(l + 0(l/n)), 

(6.45) 

ujAi{uO = ^|^e-"^W(l + 0(l/n)), -iu^Ai'{ujO = 0ie-"^W(l + 0(l/n)). 

Corresponding results can be given for Ai{uj^*), Ai'((^^*), Ai(c;j^^*) and Ai'((^^^*). Substituting (6.22), 
(6.9) and the above asymptotic formulas into (6.44), and taking into account (4.36), we readily see that 

= {-l)^-''N{z)h{z)'''-'/^e''^^-^*'>''-'h{z)-''--'/^N{z)-\l + 0(l/n)] 

= I + 0{l/n). (6.46) 



For z gTsD C+, we again have Jr{z) = I and it follows from (6.38), (6.28) and (6.31) that 

M2{z)Mi{z)-^h{zY''/^ 



X e"<^(^)'^3/i(z)-'^3/2^(^)-i_ (Q 47^ 

Moreover, we note from (6.32)-(6.33) , (5.65) and (5.48) that 

M2{z)Mi{z)-^h{zY^I'^ = e-^™Z)(z)-'^^ (6.48) 
and by (6.16) and (6.9) we obtain 

En{z) f "'"'"^f^ ~f = ^N{z)h{zrl^[I + 0(l/n)]e-^~(^)-. (6.49) 

y -u:A\{u}i) Ai(e) J 0r 

A combination of (6.47)-(6.49), (4.43) and (5.47) yields 

Js{z) = (-l)"iV(z)/l(z)'^3/2[j^Q(^/^)]g-n^(^)a3g-7V^i^a3^(^)-a3 

X e"'^(^)'^3/i(z)-'^--^/2^(z)-i 
= / + 0(l/n). (6.50) 

Let us now consider z in the interval (0, 1) on the real line. For simplicity, we only consider the case 
z £ (0, xq]. For z = X £ (0, xi], we note by (5.65) that D{z) = D*{z), where xi is the point of intersection 
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of the vertical line with the real axis; see Figure 7. It follows from (3.8)-(3.9), (6.31) and (6.38) that 



J5(x) = A^(x)e~"'^-(^)"^ 



D* (x) 



Uf=0^{x - XN,j) 



0"3 



1 

4cos^(A^7rx) 



\ciw{x) Uf=0^ix- XN,j) 



DUx) 



-as 



N{x) 



N{xy^. 



(6.51) 



4cos^(iV7r3;) n{<p++<l,-) r,* /n* M<^--<f>+) 

^w{x)D\{x)D*^{x) +' 

From (5.48), one can see that DX/D*_ = e^^'^"' and DX{x)D*_{x) = Acos^{NTTx)D{xf . Since 4){z) is 
analytic for z G (0, a), we obtain from (4.43) 



0^ 



(6.52) 



In fact, it is easily verified that D{x) is bounded and w{x)' -1 = 0{N"+f^) in this case. Since (j){x) < 
for X G (0,xi]; see (4.20) and Figure 2, e2'^^(^)tt;(3;)-iZ)(x)'2 = ©(iV^+Z^e^"^) ~ as n oo. Thus, 
Js{z) = I + 0{l/n). 

Finally, we consider the case z = x £ {xi,xo]- By (6.31), (6.33) and (3.9), we have from (6.38) 

2cos(A^7rx) ' ^ 



Jsix) = Enix)P-ix] 



{cih{x)y/^l\j=o {x-XN,j) 
2 cos(A^7rx) 



4cos^(iV7rx) 



\ciw{x)Ylf=0^{x - XN,jf 
P+{xy^En{x)-^ 



(ci/i(x))V2J]J^pl(x-3;Arj 



En{x)P^{x) 



P+{xy^En{x)~\ 



1 0^ 

h{x)w{x)^^ 1 

Since by (5.32) w{x) ~ h{x) for x G {xi,xo], it follows from (6.30) that 



(6.53) 



Jsix) = En{x)P_{x) 

= En{x)P^{x) 



1 0' 

h{x)w{x)~^ 1 



1 

-1 1 



P-{x)-^En{x)-^ 
P^{x)-^En{x)-^ 



1 

h{x)w{x)~^ — 1 1 

= /+(l/n). (6.54) 

In a similar manner, we can prove that for z in other parts of the contour F, the jump matrix Js{z) 
also tends to the identity matrix. Hence, Js{z) = I + 0(l/n) as n — )• oo on the whole contour F. By 
Theorem 3.8 in [12], the solution of the RHP for S has the asymptotic behavior 

S{z) = I + 0{l/n) (6.55) 

as n — )■ oo uniformly for z G C \ F. 
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7 Main Results 



Let Xij denote the element in the matrix X. For convenience, we define m{z) := Nii{z) and 

m*{z) := —iNi2{z)h^^{z). By (5.42) and (5.34), a straightforward calculation gives 



m{z) 



and 



m*{z) :- 



{z - a)^' + {z- byl' I Z + C/2 + {z - aY'\z - h) 
2(z-a)V4(^_ 6)1/4 I 



Nl/2 



2z 



l-z + c/2- (z-a)^/^(z-b)^/^ 
2(1 -z) 



(7.1) 



{z - 6)^2 _ _ q)1/2 + c/2 - (z - a)i/2(z - 6) V2 \ " / 1 _ ^ + c/2 + (z - a)i/2(z - 5)i/2 ~ 



2(z - ayl^z - 6)V4 



2z 



2(1 - z) 



(7.2) 



where a and 6 are given in (2.10) and c = n/N . Note that M(z) in (5.35) is analytic in C\ [a, h] and hence 
m(2;) is also analytic in C \ [a, b\. However, this is not true with the function m*{z), which has additional 
cuts (—00,0] and [l,+oo). We finally arrive at the following theorem: 

Theorem 7.1. Let I, II, III he the regions shown in Figure 6, and let I denote the Lagrange constant 
given in (4-2). With (p{z) and D{z) defined in (4-2) and (5.65), the asymptotic formula of the polynomial 
T^N,n{z) is given by 



71" AT,, 



(z) = e"'/2 J D{z)e~'"^^'hn{z) 



i + o(-) 

n 



+ S{n) 



(7.3) 



for z £ 1, where 6{n) = for z € I \ 0,± and 5{n) is exponentially small in comparison with its leading 
term for z G I Pi see Figure 1. More precisely, 5{n) = 0{N'^^^^°''^^ e^'^^^'^) , where Rei^(z) is negative 
when Rez is bounded away from the interval {a,b); see Figure 4- 
Let f{z) be defined as in (6.5). We have 



7rJv,n(^) = (-l)"^/^e"^/2 J4(^,n) 



l + 0(- 



n 



+ B{z,n) 



l + 0(-) 

n 



(7.4) 



for z G II, where 

A{z,n) = [n^'^f{z)]^/^[m{z) + m*{z)] \siu{NT:z)M{n^/^J{z))+cos{NTTz)B\{n^/^J{z))^ , 



and 



B{z, n) = [n '^/^f{z)Y^/\m{z) - m*{z)\ {sin(iV7rz)Ai'(?i2/3/(2)) + cos(iV^z)Bi'(n2/3/(^))} . 
Similarly, with f*{z) defined in (6.20), 



,n{z) = {-lfV^e^"HA*{z,n) 



l + 0(-) 
n 



+ B*{z,n) 



l + 0(-) 
n 



(7.5) 
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for z G III, where 

A*{z,n) = [n^/^f*iz)]^/^[m{z) - m*{z)] ^cos{N7rz)Bi{7i^/^ f* (z)) - sm{NTTz)Ai{n'^/^ f* (z))^ , 

and 

B*(z, n) = [n '^/^f*{z)]-^/\m{z)+m*{z)] {cos(iV7rz)Bi'(n2/3/*(z)) _ sin(iV^z)Ai'(n2/3/*(z))| . 

Proof. Since R{z) = {cie^^Y^/'^ S{z){cie''^y-^/'^R{z) by (6.37), we have 

Rii{z) = Su{z)Ru{z) + 5i2(z)i?2i(^)cie"' (7.6) 

and 

Ri2{z) = Su{z)Ri2(z) + Si2{z)R22{z)cie^' . (7.7) 

Prom (6.55), it follows that 

S^^{z) = 1 + 0(l/n) and Si2{z) = 0{l/n). (7.8) 
First, let us consider z € I. Recalling the definition of R{z) in (6.31), one obtains 

7V-1 

i?ii(z) = e^'l^e-^^^'^b{z) J{{z- XNj)-' [Sn{z)Nn{z) + Si2{z)N2i{z)] (7.9) 

j=0 

and 

N-l 

Ruiz) = cie"'/2e"'^(^)5(2)-i Hiz- xnj) [Sn{z)Nu{z) + Suiz)N22{z)] . (7.10) 

3=0 

Note by Theorem 2.1 and (3.1) that 

N-l 

ttnA^) = Yniz) = Hn{z) Y[{z- xn,j). (7.11) 

j=0 

From (3.3)-(3.4), we know that Hii{z) has different expressions in different parts of the region I. Prom 
(3.4), (7.9) and (7.8), it follows that 

N-l 

7TN,n{^) = Rniz) Hiz- XN,j) = e"""^ b{z)e-'^^^'^ m{z) [1 + 0(l/n)] (7.12) 

for z G I \ since m{z) := Nii{z). 

Recalling the notation ci = 2NTri in (3.9), we have from (3.3) and (7.8)-(7.10) 

^^LlJ zs^ie^^'^^'^ cosiNiTz) 

T^N,n{z) = Rll{z) \\{z - XN,j) - Rl2{z) 



j=o N7rw{z) llj=o {z - XN,j] 

'D(z)N^^(z)\l + 0(l/n)] T e"<^(^)^ 

D{z)w{z) 
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^1/2 1 e-''^D{z)Nu{z)[l + 0(l/n)] ^ e"'^(") 2e^^^'' cos{Nttz) ^^^^^^ ^ 0(l/n)] I (7.13) 



for z G In Q±. For simplicity, we only consider Hez < xq. It follows from (5.65) and (5.48) that 

7rN,n{z) = e"'/' {e-"<^5(z)A^ii(z)[l + 0(l/n)] T e""^^'^ D{zy^w{zy' N^iz) [1 + 0(l/n)]} . (7.14) 

Note that Re(^(z) < for z £ 1; see Figure 4. Moreover, it is readily seen from (5.47) that D{z) is 
bounded as n — ?• oo for z G I n il.±. From (5.32). we also have w{z) bounded as z — t- oo for z £ ID il.± 
and z ^ 0. Near z = 0, w{z)~^ = 0(iV"). Thus, in terms of the notation m(z) = Nii{z), 

7TN,niz) = e"'/2 | e-"<^5(z)m(z) [1 + 0(l/n)] + 0(Ar"e"i^<^<^W)} (7.15) 

for z G I n ri-t and Rez < xq. Similarly, we can show 

7rN,n{z) = e"'/2 |e-"'^z3(z)m(z)[l + 0(l/n)] + o(Ar/3e"R'=<^W)| (7.16) 

for z G In ri-t and Rez > xq. Coupling the above two formulas with (7.12) gives (7.3). 

Next, we consider the case z G II; see Figure 6. Let us first restrict z G IInC+. From (7.6) and (7.7), 
we obtain by a combination of (6.9), (6.28) and (6.31) 

Ruiz) = Sn{z)Ru{z) + 5i2(z)i22i(^)cie"' 

= (-l)"^/^e"'/^ ^^'°'^^''"^ ^ |[-zr^u;^Ai-(a;e) - €ia;Ai(a;0]iVn(z)[l + 0(l/n)] 

Uj=0 - ^N,j) I 

+ [r^^^Ai'(L^e) - ^^ujAi{uO]Ni2{z)h{z)-^[l + 0{l/n)]\ (7.17) 



and 



i2l2(^) = Suiz)Ri2iz) + 5i2(z)i?22(^)cie 
= (-l)"V^cie"^/2llS^^^ 



K-3Ai'(0 + ie4Ai(0]iVii(^)/i(^)[l + 0(l/n)] 



2 cos(iV7rz) 

+ [eUi(e) - r^Ai'(0]A^i2(^)[i + o(i/n)] 

Here, use has also been made of (7.8). Since ci = 2N-Ki, by (7.11) and (3.3) we have 

-ie^Nntz cos(Ar7rz) 



(7.18) 



N-l 



VrAr,„(z) = i?ll(z) JJ (^ - XN,j) - Rl2{z) 



j=0 



iV7ru;(z)n^0^(^; - xnj) 



see the first equality in (7.13). A combination of this with (7.17) and (7.18) gives 



vr;v,n(^) = (-l)"^/^e"'/2j 2cos(iV7rz) 



e^'^'^wiz)-^ 



-iC*^^ Ai{uj(,) - i^4ujA[{uj(,)]Nn{z)[l + 0(l/n)] 

+ [r^w^Ai'(a;e) - e^wAi(a;0]iVi2(z)/i(z)-i[l + 0(l/n)] 
[ir^Ai'(0 + i3m)]Nii{z)h{z)[l + 0(l/?i)] 



+ [eiAi(c) -r^Ai'(e)]iVi2(z)[i + o(i/n)] 



(7.19) 



39 



for z G II n C+. Recall the well-known formula of the Airy functions [9, (9.2.11)] 

Bi{z) = ±i he^'^'^^Aiioj^^z) - Ai(^ 



(7.20) 



Also, note that w{z) = h{z)[l + 0(1 /n)] as n ^ oo; see (5.32) and (5.34). This together with (7.20) gives 



7r^,„(z) = (-l)"V^e"'/2JA(z,n) 



l + 0(-) 

n 



+ B{z,n) 



l + 0( 



n 



(7.21) 



for z G IInC+. 

In a similar manner, one can see that 



TTN^^iz) = (-l)"V^e"'/2j 2cos(iV7rz) 



i^-4wAi'(w^0 + iC^uj^Ai{oj^C)] ■ Nn{z)[l + 0(l/n)] 



+ [-C-^u;Ai'{Lo^O + e^a;2Ai(a;2e)] • Ni2{z)h{zr'[l + 0(1/?!)] 



+ e 



-Nniz 



[ir 5Ai'(0 + i^4Ai(0]A^ii(^)/i(^)[l + 0(l/n)] 



+ [e^Ai(e) - r3Ai'(0]iVi2(^)[i + o(i/n)] 



(7.22) 



for z G Iln C-. Again by (7.20), we obtain the exactly same formula given in (7.21), thus proving (7.4). 
Following the same argument as given above, one can establish (7.5) for z G III. This completes the proof 
of Theorem 7.1. 

□ 

8 Special Cases 

In this section, we want to investigate the asymptotic behavior of the Hahn polynomials Qn{x] a, /3, iV — 1) 
for fixed values of x. First, we introduce the notation 



X ■.= Nz- 1/2. 



U) 



Then, it follows from (1.6) and (2.14) that 



Qn{x; a, /3, iV - 1) = Q„(iVz - 1/2; a, /?, iV - 1) 
iV"(n + a + /3 + l)„ 



(a + l)„(-Ar + i)„ 



i.2) 



Now, we derive from (7.3) asymptotic formulas for Qn{x; a, (3, N — 1) when x is a fixed number (i.e., 
z = 0{1/N)). For X > -1/2, i.e., z > 0, we obtain from (7.3), (5.65), (5.48) and (4.43) 



vr7v,n(2) ~ (-l)"2e"('/^~'^("^)) D{z) cos(A^7rz)m(z). 



(8.3) 
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Here, we have made use of the fact that e-"'?^±(^)D5^(z) = (-1)"2 cos(7V7rz)e-'^<^(^)D(z) . 
Substituting (8.1) into (5.46) gives 



n( \ (AT ^ e'''T{Nz + 1/2) 
Dlz) cos{N7rz) = r= 

Moreover, it is readily verified that 

A^"(n + Q + /3 + 1)„ ^ 



cos(A^7rz) 



gx+l/2p(^ + 1) 



27r(x + 1/2)^+1/2 



sin(7rj;). 



r(a + l)r(jV-n) 

^ ^n"+i/2r(Ar) 



as n — )• c«, and by (7.1) 

m{z) 



{a + l)ni-N+l)n 

(l + c)l/2 /c+l\"/c+l\^ (l + c)"+^+l/2 



2 • 2-i/2ci/2 V 2c 



2a+/3+l/2ga+l/2 



as z — )• 0. We now derive an exphcit formula for 1/2 — (j){z). From (4.43), (4.76) and (4.2), one has 

1/2 - 4){z) = 1/2 - 0+(z) - 7ri(l - ^z) = Re g+{z) 



(8.4) 



(8.5) 



(8.6) 



(8.7) 



for z G (0,a). Note that since x is fixed, by (8.1) z — )• as n — )■ oo. On account of (8.7) and (4.58), by 
letting n — 7- oo, we obtain 



1/2 - <A(z) 



-l-21og2 + (l + -)log(l + c) 
c 



+ z 



1 w ^ 2 1 

- iog(2;j log c 

c c c 



+ 0(z2), (8.8) 



and so 



^n{l/2-4>{z)) _ -n^-2n 



e— 2-2-(l + C)"+^(X + l/2)-+l/2g— l/2^x+l/2^-2x-l^ 



(8.9) 



(The last two equations have also been used in [10, (6.27) & (6.28)].) A combination of (8.3), (8.4), (8.6) 
and (8.9) gives 

Q„,(.; a, /3, iV - 1) ~ + ,y^N^^^,^r/.^.^n^.^^^^^ ^ ^^^^^^ 

for X > -1/2. 

Next, we consider the case when x < —1/2, i.e., z < 0. From (4.1) and (4.2), it can be shown that 
Qri(t>{z) analytic in the interval (— oo,0). Hence, we obtain from (7.3) and (5.65) 



7r;v,n(2)~I?*(^)e"('/2-<^(^»m(z). 
Note that e"('/2-<^(^)) = e"(V2-9i+(;j)) ^ ^ng+(z) ^ account of (4.58), we have 



5'+(^) 



vri - 1 -21og2+ (1 + -)log(l + c) 

c 



+ z 



- log(— zj log c 

c c c 



(8.11) 



+ 0(^2) (8.12) 



as z — )• 0. Letting n — )■ oo, we obtain 

gn(Z/2-0(^)) _ (_l)n2-2"(l + c)"+^e-"c-2^-liV-^-^/2 _ -^^2)^+ V2g-x-l/2 . 



(8.13) 
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equations (8.12) and (8.13) are also given in [10, (6.33) & (6.34)]. Inserting (8.1) in (5.49) yields 

A/2^e^+i/2 

"•W= (-.-W^./.r(-.) - 

Hence, it follows from (8.2), (8.11), (8.6) and (8.13)-(8.14) that 

Q„(.;a,/5,iV - 1) ~ ^(^^.^3^^^^ + (3.15) 

for X < -1/2. 

To conclude this chapter, we make a comparison of our results with those presented in [8] for the 
discrete Chebyshev polynomials tn{z,N). First, note that by taking a = /? = 0, the weight function 
p{x;a, P, N) in (1.2) becomes 1. Furthermore, we have 



tniz, N) = {-ir{N - n)Mz; 0, 0, iV - 1); 



.16) 



see [2, p. 174 &; p. 176] and [11]. For convenience, we consider only the case when z is real and Rez < xq- 
By (7.1) and (7.2), we have 



(z-a)V2 + (z-6)V^ {z- 6)V2 _ (^ _ a)V2 

m(z) = J J — and m [z) - 



2{z-a)i{z-b) 



2{z-a)i{z-b) 



i.l7) 



If '^N,niz) denotes the monic polynomials associated with the discrete Chebyshev polynomials, then we 
have the relationship 

1 n'^ 1 

^nANz --) = -J^UNz -^,N) = N^ttnA^), 

in view of (8.2), (8.16) and [8, (2.1)]. Since the function g{z) in (4.1) does not depend on a and /?, it is 
the same (7-function defined in [8, (3.4)]; also see (4.58) and [8, (6.25)]. For z G II (that is, xi < z < xq), 
in the notations of this paper we can rewrite (5.1) in [8] as 



7r;v,n(^) ~ (-l)"^^Fe"'/2 

X J [sin(iV7rz)Ai(n2/37(z)) + cos{NTTz)Bi{n'^/^f{z)) 



(z-a)V4 



n^/V(.)]V4 



+ 



sin(iV^z)Ai'(n2/3/(z)) +cos(iV7rz)Bi'(n2/3/(^))] |^ _ ^j!/4 [^'^'/(^)]"'^' I 



.18) 



on account of (3.37) in [8]. In view of (8.17), it can be readily seen that (8.18) agrees with (7.4). Moreover, 
we obtain from [8, (6.2)], [8, (3.34)], (5.49) and (8.17) 



„,/2^_„^(,) V2^e^%-Nzr^^ (^_a)i/2 + (,_ 5)1/2 



r(-iVz + l/2) 2(z - a)i/4(^ _ 5)1/4 
for z < 0. Similarly, by [8, (6.1)], [8, (3.33)], (5.46) and (8.17) we have 

7TN,n{z) ~ (-l)"2e"('/2-0(^))z)(^)cos(iV^z)m(z) 



£)*(^)ena/2-</'W)„,(^) (8.19) 



^.20) 
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for < z < xi. A comparison of (8.19) and (8.20) with (8.11) and (8.3) implies that our results agree 
with those in [8]. In fact, one can easily show that when q = /3 = 0, we do not need to construct the 
parametrix in region I. However, R{z) defined in (6.31) for z G II can not be extended to the region I 
when Q, /3 7^ 0. This is because the function h{z) involved in R{z) is not analytic in the neighborhoods of 
and 1. 

References 

[1] J. Baik, T. Kriecherbauer, K. T.-R. McLaughlin and P. D. Miller, Discrete Orthogonal Polynomials: 
Asymptotics and Applications, Annals of Mathematics Studies, Vol. 164. Princeton University Press, 
Princeton, NJ, 2007. 

[2] R. Beals and R. Wong, Special Functions: A Graduate Text, Cambridge University Press, Cambridge, 
2010. 

[3] D. Dai and R. Wong, Global asymptotics of Krawtchouk polynomials - a Riemann-Hilbert approach. 
Chin. Ann. Math. Ser. B 28 (2007), 1-34. 

[4] P. Deift and X. Zhou, A steepest descent method for oscillatory Riemann-Hilbert problems, Asymptotic 
for the MKdV equation, Ann. Math. 137 (1993), 295-368. 

[5] A. S. Fokas, A. R. Its and A. V. Kitaev, The isomonodromy approach to matrix models in 2D quantum 
gravity. Comm. Math. Phys. 147 (1992), 395-430. 

[6] S. Karlin and J. McGregor, The Hahn polynomials, formulas and application, Scripta Math. 26 
(1961), 33-46. 

[7] A. B. J. Kuijlaars and W. Van Assche, The asymptotic zero distribution of orthogonal polynomials 
with varying recurrence coefficients, J. Approx. Theory 99 (1999), 167-197. 

[8] Y. Lin and R. Wong, Global asymptotics of the discrete Chebyshev polynomials, preprint 

[9] F. W. J. Olver, D. W. Lozier, R. F. Boisvert and C. W. Clark, NIST Handbook of Mathematical 
Functions, Cambridge University Press, Cambridge, 2010. 

[10] Chunhua Ou and R. Wong, The Riemann-Hilbert approach to global asymptotics of discrete orthogonal 
polynomials with infinite nodes. Anal. Appl. 8 (2010), 247-286. 

[11] J. H. Pan and R. Wong, Uniform asymptotic expansions for the discrete Ghebyshev polynomials. 
Stud. Appl. Math. 128 (2012), 337-384. 

[12] W.-Y. Qiu and R. Wong, Asymptotic expansions for Riemann-Hilbert problems. Anal. Appl. 6 (2008), 
269-298. 

[13] 1. 1. Sharapudinov, Asymptotic properties of orthogonal Hahn polynomials in a discrete variable. Mat. 
Sb. 180 (1989), 1259-1277 



43 



[14] G. Szego, Orthogonal Polynomials, 4th edition, AMS Colloquium Publications, Vol. 23, Amer. Math. 
Soc, Providence R.I., 1975. 

[15] X. S. Wang and R. Wong, Global asymptotics of the Meixner polynomials, Asymptotic Analysis 75 
(2011), 211-231. 

[16] Z. Wang and R. Wong, Uniform asymptotics of the Stieltjes-Wigert polynomials via the Riemann- 
Hilbert approach, J. Math. Pures Appl. 85 (2006), 698-718. 

[17] R. Wong, Asymptotic Approximations of Integrals, Academic Press, Boston, MA, 1989. Reprinted 
by SIAM, Philadelphia, PA, 2001. 

[18] R. Wong, Lecture Notes on Applied Analysis, World Scientific, Hackensack, NJ, 2010. 



44 



I 



II ' m 



/+(()) 



